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Arctic  Mechanical  Energy  Budget 


1 .  Introduction 

The  transfer  of  mechanical  energy  from  the  atmosphere  into  the 
ocean  is  responsible  for  the  bulk  of  waves  and  currents  observed.  In 
the  Arctic  the  rate  of  transfer  will  be  greatly  altered  by  the  presence 
of  an  ice  cover.  It  is  useful  to  have  an  estimate  of  the  energy  available 
to  drive  currents  in  the  Arctic  when  constructing  large  scale  circulation 
models.  We  shall  see  that  the  ice  cover  greatly  reduces  the  available 
energy  in  some  cases.  Energy  which  is  transferred  to  the  ice  is  respon¬ 
sible  for  its  motion  and  deformation.  Both  are  important  considerations 
for  shipping,  military  operations  and  resource  recovery  in  the  Arctic. 

Noise  generation  by  ice  deformation  can  be  directly  related  to  mechanical 
energy  available. 

Mechanical  energy  flows  from  the  atmosphere  into  the  ocean.  The 
source  of  this  energy  is  solar  radiation  which  drives  the  global  atmos¬ 
pheric  heat  engine.  Work  is  done  on  the  atmosphere  and  large  scale 
pressure  gradients  develop.  The  pressure  gradient  potential  energy  may 
thus  be  taken  as  a  starting  point  in  tracing  the  flow  of  atmospheric 
mechanical  energy. 

In  the  upper  atmosphere,  far  from  the  effects  of  the  surface,  the 
balance  of  pressure  gradient  and  Coriolus  force  causes  a  geostrophic 
wind.  At  the  ocean  surface  winds  must  be  equal  to  the  current  and 
thus  small.  The  transition  region  is  a  layer  of  variable  height  known 
as  the  planetary  boundary  layer  (PBL) .  Energy  is  transferred  from  this 
layer  through  its  bottom  boundary  to  the  surface  of  the  ocean. 

Energy  may  be  transferred  to  waves  on  the  surface  or  to  ice  motions. 
Waves  may  also  be  generated  in  the  ice  cover.  In  the  treatment  that 
follows  the  word  "layer"  is  used  to  define  the  mode  of  mechanical  energy 
being  considered.  Thus,  mechanical  energy  carried  by  the  ice  is  referred 
to  as  being  contained  in  the  ice  layer  while  the  wave  mechanical  energy 
is  contained  in  the  wave  layer.  In  the  latter  case  this  is  somewhat  of 
a  misnomer  since  the  kinetic  energy  of  a  wave,  though  dropping  exponentiall 
extends  to  the  bottom.  None-the-less,  it  is  useful  to  think  of  the  ice 
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and  waves  as  surface  layers  which  mediate  the  flow  of  mechanical  energy 
into  the  ocean. 

A  schematic  drawing  of  the  energy  flow  paths  is  shown  in  Figure  1. 
Some  of  the  energy  in  the  PBL  is  transferred  directly  to  currents  in  the 
top  layer  of  the  ocean.  This  layer  has  a  strong  similarity  to  the 
atmospheric  PBL  (Coantic,  1975).  It  is  defined  by  the  depth  to  which 
the  surface  effects  are  important  and  is  marked  by  a  uniform  temperature 
and  salinity.  This  region  is  referred  to  as  the  oceanic  mixed  layer 
(OML) . 

Both  wave  energy  and  ice  M.E.  may  be  transferred  to  currents  and 
turbulence  within  the  OML.  The  energy  contained  in  the  OML  is  eventually 
transferred  to  the  rest  of  the  ocean.  This  study  seeks  to  understand 
the  effect  of  the  ice  cover  on  the  mechanical  energy  budget  of  the  Arctic. 
The  effect  will  be  limited  to  the  interface  region  from  the  top  of  the 
PBL  to  the  bottom  of  the  OML  since  these  are  the  regions  of  the  atmosphere 
and  ocean  influenced  by  the  interface. 

Each  layer  may  be  characterized  by  a  mechanical  energy  balance 
equation.  For  a  stationary  system  the  rate  of  energy  input  to  each 
layer  must  equal  the  output.  The  energy  balance  equations  are  coupled 
by  terms  describing  transfer  from  one  layer  to  another.  Other  terms 
describe  energy  transport  and  dissipation. 

Tha  balance  and  flow  of  mechanical  energy  may  be  most  conveniently 
displayed  in  a  flow  chart.  This  is  also  a  useful  way  to  represent 
the  magnitudes  and  importance  of  energy  transfer  terms  (Holland  and 
Lin,  1975).  Measurements  of  mechanical  energy  parameters  are  generally 
made  at  one  point  and  averaged  for  some  time.  Mechanical  energy  inputs 
for  ocean  models  requires  spatial  averages  on  the  model  grid  scale  and 
temporal  averages  over  the  model  time  increments.  It  is  useful  to  know 
the  geographic  variation  of  the  energy  budget  term  as  well  as  seasonal 
and  extreme  values.  All  these  considerations  lead  to  spatial  and  temporal 
averaging  scales.  These  allow  a  useful  comparison  of  mechanical  energy 
flow  in  different  regions  of  the  Arctic  at  different  times. 

The  horizontal  scale  should  be  small  enough  to  resolve  the  variety 
of  features  shown  in  Figure  2.  It  should  distinguish  between  open 


Geostrophic  Wind 


Schematic  Drawing  of  Mechanical  Energy  Transfer  Pathways  from  the 
Atmosphere  to  a  Partially  Ice  Covered  Ocean. 
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Figure  2.  Important  Features  to  be  Resolved  in  an  Arctic  Mechanical  Energy  Budget. 
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water,  pack  ice  and  the  marginal  ice  zone  (MIZ)  as  well  as  regions 
of  intense  ice  deformation.  Furthermore  winds  and  currents  should  not 
vary  appreciably  over  the  length  scale.  As  long  as  the  scale  is  less 
than  the  scale  of  synoptic  pressure  variations  this  should  be  the  case. 
Pressure  systems  in  the  Arctic  have  scales  of  500  km  or  greater  so  that 
winds  will  not  vary  much  on  a  scale  of  100  km  as  is  observed.  The  MIZ 
zone  is  about  100  m  wide  while  ridging  may  be  confined  to  a  region  a 
few  hundred  meters  wide.  If  the  scale  becomes  too  small  however,  one 
can  no  longer  see  the  forest  for  the  trees.  Ice  floes  are  10  km  or  so 
in  diameter.  Hibler  (1980)  holds  that  averaging  ice  strains  over  scales 
greater  than  20  km  averages  out  discontinuities  such  as  leads  and  shear 
ridges.  An  averaging  length  of  ~  100  km  should  thus  be  adequate  for  the 
purposes  of  this  work. 

Vertical  scales  are  determined  by  the  height  of  the  PBL  and  the 
depth  of  the  OML,  Figure  1.  These  are  of  the  order  1000  m  and  50  m 
respectively  (Coantic,  1975)  though  considerable  variation  is  seen  in 
both  layers  both  in  time  and  space.  The  vertical  scale  is  thus  much 
smaller  than  the  horizontal.  Terms  in  the  energy  balance  may  now  be 
calculated  by  integrating  the  relevant  parameters  over  the  height  of  the 
various  layers  shown  in  Figure  1.  Energy  densities  are  thus  expressed 
in  terms  of  energy  per  unit  area. 

There  are  essentially  two  seasons  in  the  Arctic,  winter  when  the 
ocean  is  almost  entirely  ice  covered  and  summer  when  large  areas  are  ice 
free.  Seasonal  variations  will  be  described  by  averages  taken  over 
periods  of  20  to  60  days.  This  is  long  enough  to  smooth  out  daily 
fluctuations  due  to  storms  and  still  short  enough  to  resolve  the  season. 
Extreme  values  due  to  storms  are  of  interest  as  well.  These  may  be 
resolved  by  taking  averages  over  a  day  since  synoptic  events  last  for 
periods  of  a  few  days. 

In  this  paper  the  mechanical  energy  budgets  for  the  PBL,  ice  pack, 
waves  and  OML  are  presented  and  integrated  over  the  appropriate  spatial 
and  temporal  scales.  Magnitude  estimates  of  the  various  terms  in  these 
budgets  are  made  based  on  the  parameters  compiled  in  Appendix  1.  These 
estimates  vary  with  geography  and  season.  A  final  section  uses  a  flow 
chart  format  to  make  comparisons  between  energy  flow  in  the  two  Arctic 
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seasons  and  for  open  water,  marginal  ice  zone  and  central  ice  pack. 

These  comparisons  give  an  insight  into  the  effect  of  sea  ice  on  mechanical 
energy  transistor  processes  and,  hopefully,  will  provide  future  workers 
with  a  qualitative  guide. 

The  main  source  of  energy  for  oceanic  circulation  is  the  atmosphere. 
Some  additional  energy  is  supplied  by  tides.  Tides  are  an  important 
energy  source  for  currents  in  shallow  waters  which  are  sheltered  from 
the  atmosphere  by  an  ice  cover.  They  are  responsible  for  tidal  cracks 
in  the  fast  ice.  However,  tidal  components  in  the  Arctic  are  small, 
(Callaway,  1976)  vanishing  at  the  pole.  For  the  most  part  they  provide 
an  insignificant  background  kinetic  energy  which  is  uniform  in  space  and 
time  for  averaging  periods  of  12  hours  or  more.  The  various  mechanisms 
by  which  tidal  energy  may  transferred  to  eddies  and  internal  waves  are 
beyond  the  scope  of  this  work. 

2.  Mechanical  Energy  Balance  in  the  PBL  and  OML 

The  mechanical  energy  balance  for  the  PBL  and  OML  may  be  derived 
from  the  classical  Navier-Stokes  force  balance  equation  for  a  fluid.  If 
the  fluid  velocity  is  u  , 

0  u  +  2pftxu  +  pVy  +  Vp  =  V*o  (2.1) 

where  2  is  the  earths  angular  velocity  which  gives  rise  to  the  Coriolus 
force,  x  is  the  gravitational  potential,  p  is  the  pressure  field  and 
o  is  the  viscous  stress  tensor 


0  =  2uD 


where 

D  =  |  (Vu  +  VCu) 


(2.2) 

(2.3) 


is  the  stretching.  Taking  the  scaler  product  of  this  equation  with 
velocity  gives  a  balance  equation  for  kinetic  and  potential  mechanical 
energy. 


(2.4) 


i  ***►■• 
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The  change  in  potential  energy  all  goes  into  gravitational  work. 


dX 

p  dl  = 


(2.5) 


the  other  terms  express  the  work  done  against  pressure  gradients 


u*7p 


(2.6) 


and  viscous  forces 


=  u*(V*a)  .  (2.7) 

In  turbulent  systems  such  as  the  atmosphere  or  ocean  it  is  necessary  to 
describe  high  frequency  velocity  components  statistically.  The  velocity 
vector  may  be  described  by  its  mean  and  fluctuating  parts 


u  =  U  +  u' 


(2.8) 


In  compressible  gases  it  is  best  to  define  the  mean  as  a  mass  weighted 
average, 

pu 

V  =  ZT  (2.9) 

P 

where  the  bar  indicates  the  conventional  averaging  operator.  The  mass 
weighted  averaging  leads  to  a  clear  separation  between  the  mean  and 
turbulent  kinetic  energies 


1  - 

2  Pu  *  u 


j  pU*U  +  j  pu’ • u * 


(2.10) 


The  mean  mechanical  energy  balance  equation  may  be  obtained  from 
the  scaler  product  of  U  and  the  mean  momentum  balance  equation, 

_  dU 

p  =  -  pVy  -  2pfbdJ  -  Vp  +  V*o  +  V* t 


(2.11) 
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which  is  essentially  the  same  as  the  total  momentum  balance  equation 
except  for  the  term  containing 


=  -  pu’  x  u'  (2.12) 

representing  turbulent  Reynold  stresses.  The  mean  mechanical  energy 
balance  equation  is  then 

C  (  4  u*u)  =  -  pU*g  -  U*Vp  +  U*(V»C)  +  U*(v*t)  .  (2.13) 

The  last  term  may  be  rewritten  using  the  chain  rule 

U*(V*T)  =  V*(t*U)  -  tr  ID  .  (2.14) 

The  first  of  these  is  the  divergence  of  a  stress  flux  and  represents  the 
at  which  turbulence  kinetic  energy  is  transferred  into  or  out  of  a  unit 
volume  bv  the  mean  flow.  The  second  is  the  rate  at  which  turbulence 
energy  is  being  produced  by  shear  in  the  mean  flow  and  is  known  as  the 
mechanical  production  term. 

It  is  reasonable  to  assume  zero  vertical  component  for  the  long  term 
mean  flow  in  both  boundary  layers.  Then 


-  pU»g  =  0 


(2.151 


which  is  to  say  the  mean  flow  on  average  does  no  gravity  work.  The  rest 
of  the  terms  may  be  evaluated  by  integrating  over  the  volume  of  each 


layer,  V 


J'  C  (y  U*U)dV  =  -  J"  U'Vpdv  +  £  U*(V'0)dv  +  J"  V*(T*U)  - 


tr  TD  dV  (2.1b) 


viscous  forces  are  important  only  in  small  scale  motions  (Phillips,  1966) 
and  may  be  neglected  in  the  mean  balance  for  both  the  atmosphere  and  ocean 
boundary  layers. 
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U*(V*c)dv  --  0 


(2.17) 


For  the  horizontal  scale  chosen,  the  mean  velocity  may  be  considered 
to  vary  only  with  height.  For  the  purposes  of  this  estimate  it  is 
assumed  that  the  turbulence  stress  tensor  is  horizontally  homogeneous. 
The  mechanical  energy  balance  equation  may  then  be  written  in  terms  of 
the  energy  flow  per  unit  horizontal  area.  For  a  PBL  of  height  h  , 


h  h  h 

/  7  57  <7  VT>dz  ■  -f  L™2  +f  {h 


•l')  -  tr  TDldz 


(2.18) 


where  k  is  a  unit  vector  normal  to  the  surface.  The  Reynolds  stress  flux 
gradient  gives  the  power  input  at  the  top  and  bottom  of  the  layer 


(k*i*U)dz  =  T,  •  U.  -  T  *U 
oz  ~  -  ~h  ~h  ~o  ~o 


(2.19) 


It  is  commonly  assumed  that  the  Reynolds  stress  at  the  top  of  the  PBL 
is  zero  (Roll,  1965).  For  the  purposes  of  estimation  it  may  also  be 
assumed  that  the  density  and  pressure  gradients  are  constant  over  the 
height  of  the  PBL  (Hunkins,  1980).  The  balance  of  mean  kinetic  energy 
for  the  PBL  may  then  be  written 


hf  <1  -T 


U)dz  =  -  /  U’Vp 


r  *U  -  I  tr  tD  dz 

-o  ~o  I 


For  a  stationary  system 


11 

37  KE  "  ■  0 


(2.20) 


(2.21) 


f  H  ,vPdz  =  -  VH0  f 


tr  tD  dz 


(2.22) 


TCS58 
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The  term  on  the  left  is  the  rate  of  working  by  the  mean  flow  on  the 
pressure  gradient.  This  is  partially  balanced  by  the  Reynolds  stress  at 
the  bottom  surface 


-  T  *U  =  +  ou'  xu'  *U 

~o  o  -  ~  ~o 


(2.23) 


where  U  is  the  velocity  of  the  water  or  ice  at  the  bottom  of  the 
~o 

boundary  layer. 

The  last  term  is  the  vertical  average  areal  density  of  rate  of  tur¬ 
bulence  kinetic  energy  production.  This  is  the  mechanical  production  term 
averaged  over  the  height  of  the  PBL.  This  turbulent  kinetic  energy  eventually 
travels  through  a  turbulence  energy  cascade  whereby  the  wavelengths  of 
disturbances  become  progressively  smaller  (Lumley  and  Panofsky,  1965). 
Eventually  the  scale  of  motion  becomes  small  enough  that  the  energy  is 
dissipated  by  viscous  heating.  Some  of  the  turbulent  kinetic  energy  may 
be  transferred  to  waves  and  trubulence  in  the  oceanic  mixed  layer  through 
work  done  by  turbulent  pressure  variations  on  the  sea  surface.  Large 
scale  pressure  variations  also  provide  an  energy  source  for  the  ocean 
but  these  are  not  associated  with  boundary  layer  transfer. 


Mechanical  Energy  Balance  in  the  OML 

The  mean  mechanical  energy  balance  for  the  mixed  layer  may  be 
handled  in  the  same  way  as  for  the  PBL.  Vertical  integration  now  extends 
from  the  bottom  of  the  layer  to  the  surface  or  bottom  of  the  ice.  The 
mean  flow  is  considered  to  be  essentially  horizontal  (Phillips,  1966) 


d_ 

dt 


(-j  pU *U)dz 


(k*fU) 


tr  tD|  dz 


(2.24) 


One  of  the  characteristics  of  the  mixed  layer  is  a  uniform  density.  The 
pressure  gradient  is  related  to  the  sea  surface  tilt  VH 


Vp  =  pgVH 


(2.25) 
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and  is  uniform  with  depth.  Assuming  no  turbulence  at  the  bottom  of 
the  mixed  layer  gives  zero  Reynolds  stress  here.  The  energy  balance 
becomes 


The  stationary  balance  is 

0  0 

U'PgVH^2  +  £  tr  xD  dz  (2.27) 

-d  -d 

The  surface  stress  power  input  is  balanced  by  changes  in  sea  surface 
tilt  and  by  turbulent  dissipation  of  energy.  In  open  water  the  surface 
stress  is  just  the  Reynolds  stress  due  to  shear  in  the  atmosphere. 

These  stresses  must  be  equal  and  opposite  in  direction  for  equilibrium 
to  occur.  Underneath  the  ice,  the  stress  is  due  to  drag  by  relative 
motion  of  the  ice  sheet. 

Sea  surface  tilt  may  be  either  a  source  or  sink  of  energy.  The 
tilt  due  to  geostrophic  pressure  variations  does  work  on  the  mean  flow. 

On  the  other  hand  some  of  the  energy  of  currents  is  transferred  tem¬ 
porarily  into  changes  in  sea  surface  potential. 

Turbulence  is  generated  in  the  OML,  as  in  the  PBL,  by  shear  in  the 
mean  flow.  Other  sources  of  turbulence  in  the  mixed  layer  are  direct 
pumping  of  turbulent  motions  by  atmospheric  turbulence  and  turbulence 
generated  by  breaking  waves.  The  first  of  these  is  thought  to  be  small. 
Breaking  of  waves  may  represent  a  large  part  of  the  turbulence  energy 
flux  to  the  ocean.  The  turbulence,  however,  is  limited  to  the  near 
surface  within  a  distance  comparable  to  the  wavelength  of  the  breaking 
wave  (Phillips,  1966).  This  dissipation  process  is  considered  in  the 
wave  energy  balance. 

3.  Waves  in  Sea  Ice 

Waves  in  open  water  near  an  ice  pack  will  transfer  some  of  their 
energy  to  the  ice  by  propagating  into  it.  Long  period  waves  may  pro¬ 
pagate  hundreds  at  kilometers  into  the  pack  (Robin,  1963).  Wadhams 
(1973)  gives  an  expression  for  the  ratio  of  transmitted  to  reflected 
wave  energy  density  for  waves  impinging  on  an  ice  sheet. 
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R 


1  + 


3  4 

32  Eh  - 

4  ? 

3og\  (l-v“) 


(3.1) 


where  h  is  the  half  thickness  of  the  ice 
E  is  Youngs  modulus 
and  \  is  the  wavelength  . 

Long  periods  are  most  effective  at  penetrating  the  ice  with  essentially 
all  the  energy  of  a  16  s  wave  being  transmitted  into  2.5  m  thick  ice. 

The  ratio  goes  as  so  that  only  25  percent  of  the  energy  of  5  s  waves 

penetrates. 

The  periods  of  waves  which  can  be  generated  by  a  storm  are  limited 
by  the  fetch  over  which  the  wind  blows.  Waves  with  periods  longer  than 
5  require  fetches  of  a  hundred  kilometers  or  more  to  gain  significant 
energy.  Open  leads  and  polvnyas  in  the  Arctic  have  fetches  which  are 
limited  to  a  few  kilometers  except  for  the  large  semipermanent  polynyas. 

For  the  most  part,  storms  over  open  leads  or  polynyas  should  not  generate 
long  swells.  The  wave  energy  generated  within  a  lead  will  thus  be 
reflected  and  scattered  by  the  ice  margin. 

It  is  possible  that  storms  over  large  open  fetches  of  water  such 
as  the  Chukchi  or  Norwegian  seas  could  generate  significant  long  period 
swells  which  would  in  turn  transfer  large  amounts  of  energy  to  the  marginal 
sea  ice  zone. 

When  a  wave  swell  travels  under  the  ice,  part  of  its  energy  is 
taken  up  as  a  flexural  gravity  wave  in  the  ice  sheet.  Attenuation  of 
this  wave  occurs  due  to  creep  in  the  ice  and  perhaps  by  fracture  as  well. 
Wadhams  (1973)  derives  the  rate  at  which  the  wave  energy  decays  with 
penetration  into  the  sheet.  Assuming  a  flow  law  for  ice  of  the  form 


dt 

dt 


n-1 


o’ 


(3.2) 


where  c  is  the  strain  tensor 

o'  is  the  deviatoric  stress  tensor 
T  =  [o'  *o'  |  is  the  effective  shear  stress 

B  =  some  function  of  temperature 
n  =  a  constant  usually  taken  to  be  3 


4 
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And  given  a  wave  profile 


f.  =  A  sin  2tt  (x/X  -  t/T)  (3.3) 

for  a  wavelength,  X  ,  and  period,  T  .  The  rate  at  which  energy  is 
lost  in  an  ice  sheet  of  thickness  h  per  unit  area  is 


.  2_ir|/,2  2.  n+1  ,n+2 

d£  _  7  4tt  ElCl/X  (1  -  v  )  h 

dt  (2B)n  (n  +  2) 

with  E  =  Young’s  modulus  for  sea  ice  -  6  GPa 
V  =  Poissons  ratio  =0.3  (Lavrov,  1969) 
h 


(3.4) 


and  B  = 


Bdz  is  the  vertical  average  of  B  .  Wadhams  estimates 


B  from  studies  of  wave  attenuation  to  be  between  107  and  10®  s-  . 

m 

The  energy  loss  may  be  averaged  over  an  entire  wave 


d£  =  ,  2tt2EA/X2  (1  -  v2)  n+1  hn+2  ,  . 

dt  "  “  (2B)n  (n  +  2) 

An  incident  wave  height  of  10  m  represents  an  upper  bound  on  the 

largest  swells  which  are  liable  to  enter  the  ice.  For  a  16s  wave,  X  =  400  m. 

0 

The  energy  dissipation  rate  at  the  ice  margin  is  500  watt/m".  For 

2 

comparison  a  1  m  wave  dissipates  only  0.5  watt/m  .  Dissipation  rates 
fall  off  rapidly  with  distance  from  the  ice  edge.  Observations  of  waves 
with  periods  from  10-60  seconds  more  than  100  km  from  the  ice  margin 
give  maximum  wave  heights  of  1  cm  for  a  dissipation  rate  of  5  x  10  ^  watt/m 

4 .  Mechanical  Energy  Balance  for  Waves 

At  the  free  surface  of  an  incompressible,  irrotational  fluid 
the  dynamic  boundary  condition  is. 


z£  = 

p  9t 


z»n 


+  ^  u«u 


sn 

z=n 


(4.1) 


Flow  Research  Report  No.  178 
December  1980 


where  C  is  a  velocity  potential, 

u  = 


(4.2) 


satisfying  Laplaces  equation. 


2 

v  0  =  0 


(4.3) 


and  r  is  the  elevation  of  the  free  surface  above  some  reference  height. 
Let  a  dot  stand  for  differentiation  with  respect  to  time, 


3n 

3t 


(4.4) 


Multiplying  Bernoulli's  equation  with  the  surface  velocity,  r,  and 
averaging  over  a  unit  area  gives  a  balance  equation  for  transfer  of 
energy  across  the  interface. 


-  p« 

n 


POn 


-  f 


u*u 


z=n 


pg 


nn 


z=n 


(4.5) 


The  term 


-  pn  =  R. 


(4.6) 


is  the  rate  at  which  the  atmosphere  does  work  on  the  waves  through 

the  product  of  pressure  on  the  free  surface  and  the  free  surface  velocity. 

The  rate  of  change  of  mean  kinetic  energy  per  unit  horizontal  area 
of  a  wave  is  given  by  Kinsman  (1965,  p.  524) 


•  •  1  2* 

p<t>n  +  |  ou  n  -  ^'F 


(4.7) 
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where 


F 


p4>udZ 


(4.8; 


is  the  horizontal  energy  flux  vector  representing  the  mean  energy  flux 
per  unit  horizontal  surface  length.  The  energy  balance  at  the  surface 
z  =  •  may  thus  be  rewritten 

t  =  pn  +  ?h*f  -  gnn  .  (4.9; 

The  energy  transferred  across  the  surface  by  cri  is  then  up  by  Kinetic 
energy  ,  T  ,  potential  energy 

V  =  gnn  ,  (4.11 


and  horizontal  energv  flux,  V  *F 

In  reality  the  mechanical  energy  budget  must  be  more  complex  than 
is  indicated  by  Equation  4.9.  The  driving  force  is  the  Reynolds  stress 
due  to  wind  shear.  This  stress  includes  both  work  done  by  small  scale 
pressure  and  shear  stresses  on  the  surface  Phillips  (1977)  suggests 
however  that  normal  stresses  provide  the  dominant  means  of  energy  trans¬ 
fer  to  waves.  Work  is  done  by  the  component  of  pressure  in  phase  with 
the  wave.  The  difficulty  comes  in  calculating  the  average  energv  flux. 

The  rate  of  energy  flux  to  the  waves  is  strongly  dependent  upon  the  wave 
structure  which  depends  in  turn  on  the  turbulence.  The  momentum  flux 
to  waves  may  be  estimated  from 

Tw  =  2.5  x  10-2  Dau2s2  (Phillips,  1977)  (4.1 

where  u  is  the  wind  at  some  reference  height,  usually  10  meters,  and 
2 

s  is  the  mean  square  slope  of  the  waves,  at  saturation 

s2  ~  1.4  x  10 


(4 .  i ; 
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This  gives 

T  -  3.5  x  10"4  p  u2  (4.13) 

w  —  a 

equivalent  to  a  wave  drag  coefficient 

C  =  3.5  X  10"4  .  (4.14) 

w 

This  may  be  compared  to  the  drag  coefficient  giving  the  surface  shear 
stress 


T 

O 


C  p  u 
z  a  z 


with 


-10 


=  1.5 


10 


-3 


(4.15) 


(4.16) 


The  momentum  flux  to  waves  is  thus  about  one  fifth  the  flux  to  currents. 
Mechanical  energy  transfer  may  be  much  larger  however  since  power  density 
is  the  product  of  stress  and  velocity.  In  the  case  of  waves  the  appro¬ 
priate  velocity  is  the  phase  velocity  which  is  much  larger  than  the  water 
surface  velocity. 

The  wave  energy,  E  =  T  +  V  ,  is  dissipated  by  viscosity  and  wave 
breaking.  Phillips  (1977)  gives  an  attenuation  coefficient 

3F 

Yv  =  H  2E  =  2v/k  (4.17) 


relating  the  energy  loss  rate  to  the  energy  contained  by  a  given  wave 
number,  k  ,  and  the  viscosity. 

In  order  to  evaluate  the  total  energy  dissipation  rate  integrated 
over  all  wave  numbers  the  wave  spectrum  must  be  known  as  well  as  the 
amplitude  of  the  waves.  For  the  most  part  however,  the  viscous  attenua¬ 
tion  term  is  small  especially  at  low  wavenumbers  where  the  mechanical  energy 
of  ocean  waves  is  concentrated.  Phillips  (1977,  p.  54)  gives  attenuation 


I 


-17- 
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coefficients  as  a  function  of  wave  number.  On  a  clean  surface  a  wave 

with  wavelength  X  =  1  m  or  more  has  an  attenuation  coefficient  of  less 

-4  -1 

than  2  x  10  s  .  For  this  reason  attenuation  is  not  considered 
important  in  most  descriptions  of  gravity  waves.  On  the  other  hand  it 
is  possible  for  low  frequency  waves  to  feed  energy  into  higher  frequencies 
where  dissipation  occurs. 

Other  energy  may  be  dissipated  by  wave  breaking.  It  is  difficult 
to  assess  the  magnitude  of  this  dissipation  mechanism.  It  is  most 
important  in  the  development  of  a  saturated  sea  state.  In  this  case 
wave  breaking  is  wide-spread  and  is  responsible  for  the  form  of  the  wave 
spectrum.  Small  wavelengths  break  at  amplitudes  smaller  than  large 
wave  lengths  giving  a  saturation  spectrum  which  decays  at  short  wave¬ 
lengths  (Phillips,  1966,  p.  113). 

Waves  also  transfer  energy  to  current  through  the  work  of  the 
radiation  stress  tensor.  Its  magnitude  is  related  to  the  energy  density 
of  the  wave  bv 


S11  =  \  E(K)  (Kraus,  1972) 


(4.18) 


The  rate  at  which  the  radiation  stress  does  work  on  a  current  is  the 
product  of  radiation  stress  and  current  velocity  gradient.  The  energy 
density  of  the  waves  also  increases  as  they  pass  into  an  opposing  current 
by  geometric  covergence.  A  wave  of  velocity  Cq  and  energy  density 
which  moves  into  a  current,  u,  has  a  new  phase  speed 


and  energy 


C  = 


C 

o 


(1 


o 


E  = 


E  C2 
o  o 


C(C  +  2u) 


(i.19) 


(4.20) 


The  rate  at  which  the  waves  work  on  a  current  is 


1  i  h 

(u  +  ^  C)  +  2  E  SJ  (Phillips,  1966)  (4.21) 


-3E  _  3E 
3t  3x 
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This  last  term  may  be  added  to  the  steady  state  mechanical  energy 

•  • 

balance  equation,  T  =  o  ,  V  =  o  , 

-I  n  =  V-F  +  -  (u  +  i  C)  ^  E  — £  (A. 22) 

w  H  ~  2  dx  2  dx 

note  that  P  has  been  replaced  by  as  a  more  complete  surface  stress 

representation.  Dissipation  by  wave  breaking  is  included  as  an  unknown 
function,  £  . 

5 .  Mechanical  Energy  Balance  for  Sea  Ice 

The  momentum  balance  equation  for  an  area  of  sea  ice  may  be  used 
to  derive  the  mechanical  energy  balance.  In  the  plane  of  the  sea  ice 
the  two  dimensional  momentum  balance  is  (Coon  and  Pritchard,  1979) 

mu  =  -2  m  E  A  u  +  7  •  o  +  t  +T  -  mg  7H  (5.1) 

where  m  is  the  mass  per  unit  area  of  sea  ice,  u  is  the  horizontal 
velocity  vector,  il  is  the  angular  velocity  vector  at  the  earth,  c 
is  the  Cauchy  stress  resultant  in  excess  of  hydrostatic  equilibrium 
T  is  the  traction  exerted  on  the  upper  surface  of  the  ice  by  the 
atmosphere,  is  the  traction  on  the  lower  surface  due  to  water  drag, 

g  is  gravitational  acceleration  and  H  is  the  dynamic  sea  surface 
height . 

Taking  the  inner  product  of  the  momentum  balance  with  ice  velocity 
gives  the  mechanical  energy  balance  equation 

d  1 

(y  mu*u)  =  u  •  (V*o)  +  u  •T^+u  •  T^-mgu  •  VH  (5.2) 

The  Coriolus  term  drops  out  since  the  force  is  perpendicular  to  the  ice 
motion  and  thus  does  no  work.  The  ice  is  driven  by  the  atmosperic 
power  term, 


P 


a 


u  ♦  T 

-  „a 


(5.3) 
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t 


t 


The  ice  transfers  energy  to  the  ocean  at  a  rate 

p  =  -  u  •  T  (5.4) 

W  -v  ~w 

both  by  exerting  drag  on  the  mixed  layer.  The  potential  energy  of  the 
ice  is  increased  by  work  done  against  the  dynamic  topography 

PH  =  mg  u  •  VH  (5.5) 

The  remaining  term  represents  the  energy  input  to  the  ice  itself  due 
to  divergence  of  the  stress 


I 


PT  *  -  u  •  (V  •  0) 


(5.6) 


This  term  may  be  written 


I 


♦ 


u  •  V»o  =  T  •  (u  a)  +  tr  o  D 


(5.7) 


The  first  term  is  the  divergence  of  stress  flux;  its  physical  meaning 
will  become  more  apparent  when  integrated  over  the  area  of  the  grid. 
The  second  term  is  the  trace  of  the  product  of  stress  and  stretching. 
It  represents  the  work  done  on  the  ice  due  to  both  elastic  and  plastic 
deformations.  The  stretching  D  may  be  determined  from  velocity 
gradients  within  the  grid. 


D  = 


(Vu  +  VTu) 
2 


(5.8) 


where  T  signifies  the  transpose  operation.  The  rate  of  work  on  the 
ice  is  thus. 


p.  =  tr  ~  D 

i  ~  ~ 


(5.9) 


X 
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and  the  mechanical  energy  balance  may  be  written 


IT  (2  m 


=  V.(ua)  +  p.  +  pu  +  Pw  +  PH 


(5.10) 


The  mechanical  energy  balance  for  the  ice  layer  may  be  integrated 
over  the  area  of  the  ice  within  a  grid  element.  If  we  denote  this 
a  real  average  by 


■  hf* 


the  integrated  equation  may  be  written 


3  ,1  .  1  f 

3F  (2  m  ^  =  AT  J 


V*(ua)  dA.  +  p,  +  p  +  p  +  p 
i  l  w  ra  H 


(5.11) 


(5.12) 


The  integral  of  stress  flux  divergence  may  be  evaluated  using  the 


Green-Gauss  theorem 


dA^  =  ~  £  (ua)  *n 


(5.13) 


This  term  is  thus  the  rate  at  which  work  is  transferred  to  the  area  A. 

i 

by  tractions,  c»n  ,  on  its  boundary,  F  ;  n  is  a  unit  normal  to  the 
curve  T  . 

For  the  length  scale  of  100  km  chosen  the  mean  10  m  wind  may  be 
considered  to  be  a  constant  over  the  entire  averaging  area.  This  gives 
a  constant  wind  stress  of 


Iw  =  C10  ~10 '  HlO 


(5.14) 


where  is  the  drag  coefficient  referenced  to  the  mean  10  m  wind 

velocity  .  The  mean  currents  under  the  ice  may  also  be  considered 

to  be  constant  over  the  area  of  a  single  grid. 


T  =  c  u  u 

~w  w  1  ~w 1  ~w 


(5.13) 
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The  large  scale  pressure  variations  which  give  rise  to  the  dynamic 
topography  are  large  enough  that  VH  may  be  considered  to  be  constant 
over  the  area  of  the  grid  (Newton,  1975). 

6 .  Flow  Chart  Description 

The  energy  transfer  mechanisms  associated  with  the  power  input  from 
the  atmosphere  to  the  ocean  have  been  described  by  four  mechanical 
energy  balance  equations.  Table  1.  Each  of  these  equations  is  coupled 
to  the  others;  a  power  output  from  one  equation  being  an  input  to 
another.  Keeping  track  of  four  energy  balance  equations  simultaneously 
is  unwieldy.  Mechanical  energy  tranfer  pathways  are  best  described 
by  a  flow  chart.  Figure  3. 

This  chart  represents  a  steady  state  flow  of  energy  from  the 
atmosphere  into  the  ocean.  The  energy  density  within  any  given  layer  is 
assumed  to  not  change  very  much  over  a  period  of  one  day.  Thus  the 
energy  density  may  be  estimated  by  using  daily  averages  for  the  relevant 
parameters.  Since  average  energy  density  is  essentially  constant,  for  a 
day,  its  time  derivative  is  zero.  The  remaining  terms  in  the  energy 
budget  must  then  balance  each  other.  In  other  words  the  power  input 
into  each  layer  must  equal  the  power  out.  The  flow  chart  shows  the 
various  transfer  (power)  terms  considered  here  as  inputs  and  outputs  to 
each  layer. 

Power  input  to  a  layer  can  arise  from  one  of  two  sources.  An 
adjacent  layer  may  transfer  part  of  its  energy  to  it.  Thus  the  atmosphere 
does  work  on  ice,  waves  and  currents  and  transfers  a  small  fraction  of 
its  energy  to  them.  The  transfer  is  uniformly  downwards  because  of  the 
density  and  velocity  contrasts  between  air  and  water.  Similarly  currents 
may  exchange  energy  with  the  ice  and  waves.  The  energy  flow  here  may 
change  direction,  however,  and  this  is  indicated  by  a  double  arrow  on 
the  flow  chart.  In  general  the  net  flux  of  energy  will  be  into  the 
ocean. 

A  layer  may  also  be  subject  to  an  energy  flux  due  to  transmission 
(as  opposed  to  transfer)  of  energy  within  the  layer  from  outside  the 
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Table  1.  Mechanical  Energy  Budgets  (Steady  State), 


Planetary  Boundary  Laver: 


0  =  /  L:  dz  •  u1  —I  *U  I  tr  I  D  dz 

I  ~o  ~o  I  ~  ~a 


h 

-/ 


Oceanic  Mixed  Laver: 


0  -  -  I  •  U  + 

-o  -o 


o 

/ 


pgVHdz  +  I  tr  T  dz 


Waves : 


o  ■  -  V'  vh  •  r  +  « +  i «  f  +  i  E  s 


=  ^  J' (u  c)  *n 


dl  +  u  •  T  +  u * t  +  tr  aD  +  mg  u-VH 

-  -a  ~  -v  ~~  - 


in  each  equation  the  velocities  are  given  by  u  and  tractions  at  the 
interfaces  bv  t  . 


PLANETARY  BOUNDARY  LAYER 


Figure  3.  Chart  of  Mechanical  Energy  flow  from  the  Atmosphere  into  a  Partially  Ice 
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grid.  This  is  particularly  important  in  the  ice  where  boundary 
tractions  can  cause  significant  energy  to  be  transferred  over  hundreds 
of  kilometers.  Long  period  waves  will  also  transmit  energy  over  longer 
distances  under  some  conditions.  This  energy  input  mechanism  is 
referred  to  as  lateral  flux  in  the  flow  chart.  This  flux  may  be  either 
positive  or  negative  as  indicated  by  the  double  arrows  in  Figure  3. 

No  lateral  flux  terms  are  included  for  he  two  boundary  layers.  In  both 
the  PBL  and  OML  velocity  gradients  are  much  larger  in  the  vertical 
direction  than  horizontally.  Furthermore  the  horizontal  surfaces 
through  which  energy  is  transferred  are  much  larger  than  the  lateral 
layer  boundaries.  Both  of  these  considerations  imply  that  lateral 
energy  flux  within  these  layers  is  negligible  compared  to  the  vertical 
fluxes. 

The  flow  chart  includes  a  term  for  transfer  of  energy  between  waves 
and  ice.  The  magnitude  of  this  transfer  has  been  shown  to  be  small 
with  most  of  the  wave  energy  reflected  from  the  ice.  However,  the 
fracturing  of  ice  by  ocean  swells  is  important  for  noise  generation  in 
the  marginal  ice  zone  and  for  inducing  breakup  of  the  ice  pack  in  spring. 
The  term  is  thus  included  for  comparison  at  different  localities  and 
seasons . 

Each  layer  may  lose  energy  by  transfer  or  transmission.  In  fact 
most  of  the  energy  lost  will  be  seen  to  be  due  to  various  dissipation 
mechanisms.  In  the  boundary  layers  and  waves  dissipation  occurs  when 
the  turbulence  energy  cascade  transfers  mechanical  energy  to  turbulent 
motions  of  constantly  decreasing  wavelength.  Eventually  the  wavelengths 
approach  molecular  scales  and  the  viscosity  of  the  medium  dissipates  the 
energy  as  heat. 

Deformation  of  the  ice  pack  also  dissipates  energy;  as  heat  due 
to  the  friction  of  the  ice  as  it  grinds  against  itself  and  as  potential 
energy  stored  in  ridge-keel  systems.  This  potential  energy  is  lost 
when  the  ice  melts  and  in  fact  adds  a  miniscule  amount  of  heat  during 
the  melting  process.  All  of  the  terms  considered  in  the  mechanical 
energy  budget  re  orders  of  magnitude  smaller  than  the  fluxes  of  heat 
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,  due  to  radiation  and  phase  changes  considered  in  a  thermal  energy 

budget . 

The  dissipation  terms  are  the  most  difficult  to  measure.  An 
absolute  measure  would  require  the  detection  of  small  changes  of 
^  heat  flux  in  a  system  overwhelmed  by  much  larger  changes.  Some 

measures  of  turbulence  dissipation  have  been  made  in  the  atmosphere. 

For  the  most  part  the  dissipation  term  is  considered  as  a  residual  in 
energy  balance  equations  where  all  the  other  terms  are  measured. 

I  Measurements  of  energy  dissipation  by  turbulence  in  the  upper  layer 

of  the  ocean  and  by  wave  generated  turbulence  are  essentially  non-existant 
due  to  the  limitations  of  transducers.  Again  these  quantities  are 
given  as  residuals.  It  is  possible  to  measure  the  potential  energy  of 
p  ridges  and  keels  directly  and  these  estimates  are  included  here. 

The  plastic  properties  of  ice  as  measured  in  the  laboratory  may  be  used 
in  combination  with  ice  dynamics  results  to  estimate  the  dissipation 
rates  due  to  ice  strains. 

i 


» 


» 


» 


» 
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7.  Magnitudes  of  Transfer  Terms 

The  source  of  energy  for  winds  in  the  PBL  is  the  pressure  gradient. 
This  term  gives  rise  to  a  geostrophic  wind  which  is  independent  of  the 
surface  conditions,  be  it  ice  covered  or  open  water.  Furthermore,  the 
geostrophic  wind  is  essentially  independent  of  season  in  the  Arctic 
(Albright,  1980).  The  energy  input  to  the  mean  flow  in  the  boundary 
layer  is  given  by  the  scaler  product  of  pressure  gradient  and  mean  wind 
integrated  over  the  thickness  of  the  layer. 


PG,  PBL 


(7.1) 


At  the  top  of  the  PBL  the  wind  is  geostrophic  and  perpendicular  to  the 
pressure  gradient  so  that  no  work  is  done.  The  wind  turns  as  it  approaches 
the  bottom  of  the  layer  by  an  angle,  y  ,  at  about  20°  (Carsey,  1980)  in 
the  Arctic.  The  turning  angle  is  a  function  of  the  stability  at  the 
atmosphere  varying  from  10°  to  50°  with  increasing  stability  (Brown, 

1980).  In  general  the  atmosphere  over  sea  ice  is  stable  since  the  ice 
acts  to  insulate  the  atmosphere  from  heating  by  the  ocean.  Strong 
instability  may  occur  in  winter  over  open  leads  and  polynyas  where  the 
ocean  provices  large  amounts  of  heat  to  the  atmosphere. 

The  Eckman-Taylor  model  may  be  used  to  provide  a  crude  approxi¬ 
mation  to  the  energy  delivered  to  the  PBL.  If  the  geostrophic  wind  is 
G  the  velocity  parallel  to  the  pressure  gradient  is 


u  (z)  =  G  (1-  e  KZ  cos  kz)  (7.2) 

P 

l  _ 

where  v  =  (f/2v)‘  is  the  Eckman  depth  which  depends  on  the  Coriolus 
parameter  f  and  the  eddy  viscosity  v  .  This  expression  may  be  inte¬ 
grated  through  the  Eckman  depth  for  a  reasonable  approximation  to  the 
average  velocity  parallel  to  pressure  gradient  in  the  boundary  layer. 


u 

P 


(z) 


K 


1/k 


u(z)  dz 


(7.3a) 
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G  (1-  e  cos  kz)  dz 


(7.3b) 


(sin  1  -  cos  1) 


(7.3c) 


.45  G 


(7.3d) 


The  annual  mean  geostrophic  wind  during  AIDEX  (1975-1976)  was 


G  =  8.8  m/s 


(7.4) 


which  gives 


Up (2 )  =  4.0  m/s 


(7.5) 


The  geostrophic  wind  magnitude  is  related  to  horizontal  pressure  gradient 


?.P 

p  f  !  h  1 

a 


(7.6) 


f  =  2f.  sin  if  *  -  1.4  x  10  ^  s  ^  at  76°  N  latitude 


1  3 

i.J  3 

m 


|Vh  P|  ~  1.6  x  10  Pa/m 


(7.7) 


Finally  the  work  done  on  the  PBL  per  unit  area  may  be  estimated 


G,PBL 


u(z)  dz  *Vp  ~  ku^  |Vhp| 


(7.8a) 


-6.3  k 


(7.8b) 
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typically  the  height 


-28- 


of  the  PBL  is  10' 


m  so  the  total  power  input 


is 


C,  PBL  - 


~  6 


watt 


(7.9) 


Stability  conditions  can  effect  this  value  in  two  ways.  At  low  stabilities 
the  turning  angle  is  low  and  the  pressure  gradient  cannot  do  much  work. 

On  the  other  hand,  the  height  of  the  PBL  increases  under  these  condi¬ 
tions  so  that  the  total  work  done  over  the  layer  thickness  may  be 
larger. 


Reynolds  Stress  Work 

The  work  done  by  the  PBL  on  the  surface  of  the  ocean  is  simply 
the  product  of  Reynolds  stress  and  the  appropriate  surface  velocity. 

Ice  and  current  velocities  are  typically  about  2%  of  the  10  meter  wind 
and  are  turned  at  about  40°  to  the  right  of  it  in  summer  (McPhee  1980) . 
The  velocity  component  parallel  to  the  geostrophic  wind  is  thus 


Up  -  u^0  (.02)  cos  40° 

or 

u  .015  u,~ 

p  10 

The  power  is  obtained  by  multiplying  by  the  Reynolds  stress. 


T 

a 


'10 


2 

U10 


(7.10) 

(7.10a) 


(7.11) 


=  .015 


Pa  C10  U10 


(7.12) 


Typically  u 


10 


T  u  -  5  m  watt/m 
a  p 


is  0.5  m/s  which  gives  a  transfer  power  of 
2 

This  is  just  a  small  fraction  (-.1%)  of  the 


energy  input  into  the  PBL  .  The  remainder  of  the  energy  must  be  dissipated 
by  the  turbulent  energy  cascade  within  the  PBL.  Extreme  winds  of  20  m/s 
are  occasionally  observed  in  the  Arctic.  The  power  input  to  the  ice 


K 
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input  to  the  ice  by  such  a  wind  would  be, 


T 

a 


u 

P 


300 


mwatt 

2 

m 


(7.13) 


The  energy  available  in  the  PBL  will  also  be  higher,  assuming  that 
the  geostrophic  wind  is  ~40  m/s  gives 


P„  -  25  — r—  (~1%  available  power)  (7.14) 

(j,  rDL  Z 

m 

In  the  winter  the  ice  velocity  is  about  .0065  u^q  (McPhee,  1980). 

The  typical  power  input  to  the  ice  thus  drops  to  ~2  - - —  with 

an  extreme  value  of  ~100  ■  .  When  the  ice  is  mo?ionless  of 

course  no  energy  is  transferred  to  it  and  all  of  the  energy  in  the 
PBL  is  dissipated. 

The  energy  transferred  to  currents  will  be  very  similar  to  the 
values  for  the  ice  in  summer.  The  drag  coefficient  over  water  varies 

_3 

with  the  sea  surface  state.  Kraus  (1972)  gives  C^q  -  1 . 3  x  10 
which  means  that  the  power  transfer  coefficients  are  -.65  times  the 
values  given  for  ice  in  the  summer.  The  power  input  to  the  currents 
per  unit  area  will  not  vary  much  with  time  of  year.  The  total  power 
within  a  unit  grid  available  to  drive  currents  will  be  affected  most 
by  the  amount  of  open  water  available.  In  open  water,  of  course, 
all  the  power  goes  into  the  water.  In  the  Arctic  summertime  there 
is  typically  only  10%  open  water  overall,  but  large  open  areas  of  water 
appear  occasionally. 

In  winter  much  less  open  water  is  available  to  transfer  energy 
from  winds  directly  to  currents.  Furthermore  ice  motions  are  restricted 
in  the  winter  so  that  the  ice  transfers  less  energy  to  currents.  In 
summer  free  drift  ice  motions  should  transfer  a  considerable  part  of 
their  energy  to  the  oceanic  mixed  layer.  The  bulk  of  this  energy  will 
be  dissipated  in  the  OML  as  turbulence,  but  some  will  be  transferred 
to  currents  by  the  work  of  Reynolds  stresses.  When  the  ice  moves  with 
a  current  no  work  is  done  since  the  Reynolds  stress  is  zero.  When  the 
current  is  zero  or  at  9(j°  to  the  ice  motion  no  work  is  done  either  and 
all  the  kinetic  energy  is  lost  to  turbulence.  Energy  transer  occurs 


» 


) 


» 


» 


» 


» 


» 


‘  I 


Flow  Research  Report  No.  178 
December  1980 

-30- 

when  the  ice  and  current  have  velocity  components  which  are  parallel 
but  different.  Energy  may  be  transferred  from  ice  to  currents  and 
vice  versa. 

Typical  ice  velocities  with  respect  to  currents  are  1  to  14  cm/s. 

A  lm  water  drag  coefficient  has  been  estimated  by  Langleben  (1980)  as 

_3 

=  4.1  x  10  .  The  shear  stress  is 


T 

W 


(7.15) 


where  u,  is  the  relative  current  velocitv  1  m  beneath  the  ice.  The 
1 

energy  transferred  into  or  out  of  the  ice  will  be 


or 


P. 

i,c 


i  T 

w 


(7.16) 


P. 

i,c 


1  P. 


w 


(7.17) 


u  is  typically  between  1  and  10  cm/s  giving  a  transfer  power  of  0.004 

1  2 
to  4  mwatt/m 

2 

The  energy  input  to  the  ice  has  been  estimated  at  5  to  300  mwatt/r.~ 

The  difference  between  this  value  and  the  power  transferred  to  the  0>IL 

must  be  taken  up  by  ice  deformation.  Pritchard  (1981)  has  estimated  the 

2 

rates  of  energy  dissipation  by  ice  deformation  at  0  to  200  mwatt/m*’  by 
using  the  full  AIDJEX  model.  Almost  all  of  this  power  could  be  transferred 
from  one  area  to  another  via  the  stress  flux.  The  rate  of  working  bv  winds 
upon  a  fully  developed  sea  is,  (Kraus,  1972) 


W  ~  .05  u^q  i  (7 . 18) 

2  -3 

where  T  is  the  Reynolds  stress  t  =  p  C,  „  u, n  with  C, _  =  1.3  xlO 

a  10  10  10 

The  rate  of  working  is  a  maximum  in  a  fully  developed  sea  since 
wave  velocities  and  surface  roughness  are  maximized.  In  the  Arctic  the 
fetch  is  rarely  large  enough  to  generate  such  a  sea  state  and  in  fact 
fully  developed  seas  are  found  only  occasionally  in  the  open  ocean.  The 


I 
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rate  of  working  on  intermediate  sea  states  is  not  known  or  even  esti¬ 
mated  but  must  be  less  than  the  maximum  value  given.  The  maximum  for 
a  typical  wind  of  5  m/s  is  about 

W  ~  10  — (7.19) 

max  2 

m 

If  the  wind  reaches  an  extreme  value  of  20  m/s  this  power  can  increase 
drastically  to 


W 

max 


700 


m  watt 
2 


(7.20) 


m 

Arbitrarily  assuming  half  the  power  for  an  intermediate  sea  state  gives 
rates  of  working  on  waves  roughly  equivalent  to  the  rates  on  ice  or 
currents. 


8.  Discussion 

Mechanical  energy  flow  charts  for  a  variety  of  Arctic  conditions 
are  shown  in  Figure  A.  Four  geographic  locations  are  shown.  For  two 
of  these,  Central  Arctic  Basin  (CAB)  and  Shear  Zone,  conditions  are 
significantly  different  in  summer  and  winter  and  these  are  shown. 

The  flow  of  mechanical  energy  into  the  open  ocean  is  considered 
first  as  a  reference  point.  In  the  Arctic  summer  and  winter  atmospheric 
energy  levels  are  not  significantly  different,  so  only  one  flow  chart 
is  shown.  The  mechanical  energy  of  the  PBL  is  mostly  dissipated  by 
turbulence.  The  remainder  is  shown  to  be  split  evenly  between  waves 
and  the  OML  .  These  power  levels  are  somewhat  arbitrary  since  it  is 
difficult  to  separate  the  direct  driving  of  the  OML  from  that  mediated 
by  waves.  Much  of  the  wave  energy  is  lost  due  to  wave  breaking.  Some  is 
transferred  to  the  OML  currents  and  a  considerable  amount  can  be  radiated 
out  of  the  area.  This  is  the  major  source  for  long  wavelength  wave  energy 
in  the  marginal  ice  zone  (MIZ).  Only  in  open  water  is  the  fetch 
sufficient  to  generate  significant  wave  energy. 

The  mechanical  energy  flow  in  the  marginal  ice  zone  is  considerably 
more  complex.  Assuming  50  %  ice  in  this  region  means  that  about  half  the 
available  energy  from  the  PBL  drives  the  ice  directly.  Ice  mechanical 
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energv  can  be  dissipated  by  floe  bumping  and  grinding  and  by  Reynolds 
stresses  in  the  OML  .  The  interaction  of  swells  with  the  ice  can  also 
be  significant  near  the  ice  margin.  The  source  of  wave  swell  energy 
in  the  MIZ  is  swells  generated  outside  the  area.  Short  wave  lengths 
may  be  generated  within  the  MIZ.  Their  interaction  with  the  OML  will 
be  the  same  as  in  open  water.  The  presence  of  ice  will  severely  limit 
wave  generation  by  dissipating  wave  energy  in  breaking  against  ice 
flows . 

The  ice  completely  dominates  energy  flow  in  the  CAB.  Only  1-2% 
open  water  exists  here  in  winter,  so  all  the  mechanical  energy  available 
from  the  PBL  moves  into  the  ice.  The  ice  moves  more  slowly  than  open 
water  currents  so  the  energy  transferred  is  low.  Furthermore,  ice 
deformation  tends  to  concentrate  near  shorelines  so  energy  dissipation 
should  be  low.  The  energy  transferred  to  the  ice  will  be  transmitted 
by  stress  flux  to  other  regions  especially  the  Shear  Zone.  A  small 
amount  of  energy  will  be  transferred  to  the  OML  by  the  ice  motions. 

In  the  summer  as  much  as  10%  open  water  is  found  in  the  CAB. 

Energy  transfers  power  to  waves  and  OML  currents  are  thus  shown  to  be 
about  10%  of  that  input  to  the  ice.  The  open  water  also  reduces  the 
strength  of  the  ice  to  zero  thus  essentially  no  mechanical  energy  is 
transmitted  to  the  Shear  Zone  in  summer.  The  ice  moves  more  in  summer 
so  more  energy  is  transferred  to  it.  A  considerable  amount  of  this 
energy  is  lost  by  deformation  at  flow  boundaries  in  the  summer.  There¬ 
fore,  less  energy  is  left  to  transfer  into  the  OML  than  in  open  water. 

This  is  why  Arctic  Ocean  currents  are  considerably  less  energetic  than 
their  mid-latitude  counterparts. 

Only  one  flow  chart  is  shown  for  the  Shear  Zone  in  winter.  In  the 
summer  the  mechanical  energy  flow  is  essentially  the  same  as  in  the  CAB. 

In  winter  the  only  difference  is  the  mechanical  transmitted  into  the 
Shear  Zone  from  the  CAB.  This  causes  a  large  increase  in  the  dissipation 
of  energy  by  ridging  and  other  forms  of  deformation  in  this  region. 
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9.  Conclusion 

Mechanical  energy  balance  equations  have  been  written  for  the 
various  "layers"  important  to  the  transfer  of  energy  from  the  atmosphere 
into  an  ice  covered  ocean.  These  equations  are  coupled  by  the  transfer 
of  energy  from  one  into  another.  This  coupling  is  presented  graphically 
in  the  form  of  a  mechanical  energy  flow  chart.  The  magnitudes  of  the 
transfer  terms  have  been  estimated  and  presented  in  flow  chart  form  for 
a  variety  of  Arctic  conditions. 

The  presence  of  a  sea  ice  cover  has  been  shown  to  have  a  profound 
effect  on  the  flow  of  mechanical  energy  into  the  Arctic  Ocean.  In  both 
summer  and  winter  the  energy  available  to  drive  currents  is  severely 
restricted.  Instead,  this  energy  goes  into  deforming  the  ice  pack.  Most 
of  the  deformation  is  concentrated  in  a  Shear  Zone.  The  levels  of  power 
being  dissipated  have  been  estimated. 

The  interaction  of  waves  with  ice  has  also  been  considered.  The 
presence  of  ice  provides  an  important  means  for  dissipating  wave  energy 
by  wave  breaking.  Long  period  swells  generated  in  the  open  ocean 
dissipate  their  energy  quickly  at  the  ice  margin.  The  levels  of  power 
transfer  indicate  that  this  may  be  an  important  noise  generation  region. 
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Appendix 

Typical  and  extreme  values  for  currents,  ice  motions,  winds,  and 
ice  cover  in  the  Beaufort  Sea. 

A)  Currents 

Mixed  layer  currents  to  ~  30  m  .  Mean  surface  drift  velocities 

averaged  over  20  days  during  AIDJEX  Pilot  project  were  2  cm/s  to  the 

southwest  over  the  CAB  (Hunkins,  1974). 

Currents  in  the  off-shore  region  may  be  considerably  higher  typically 

5  cm/ s  for  the  prevailing  currents  to  50  cm/s  around  points  or  driving 

storms  in  open  water  (Brower,  et  al.,  1977). 

Fluctuating  parts  of  currents  are  larger  than  the  mean  typically  10 

cm/s  in  CAB  with  extreme  to  14  cm/s  under  ice  turbulent  kinetic  energy 

3 

at  7  ergs/cm  (Hunkins,  1974).  McPhee  (1980)  equates  free  drift  of  ice 
with  mixed  layer  current  in  summer  ice  drifts  at  0-25  cm/s.  Wiseman 
(1974)  has  observed  surface  currents  averaging  14  cm/s  with  extremes  to 
34  cm/s  in  open  water,  offshore,  driven  by  winds.  Surface  currents 
under  fast  ice  in  shallow  water  are  very  small  5  cm/s  (Anguard,  1980). 

Sub-surface  flow  ~  50  to  300  m  depth.  Mean  flow  ~  same  as  surface 
drift  maybe  a  little  less  ~  2  cm/s  (Hunkins,  1974).  Fluctuations  flow  - 
baroclinic  eddies  velocity  maximum  at  ~  150  m  at  40  cm/s  or  more  is 
often  observed.  Kinetic  energy  at  this  level  is  63  erg/cm'*  (Hunkins, 

1974,  1980). 

B)  Winds 

Geostrophic  calculations  gives  synoptic  free  stream  flow  above  the 
planetary  boundary  layer  averages  ~  9  m/s  both  summer  and  winter. 

Surface  winds  due  to  geostrophic  winds  average  ~  5  m/s  summer  and  winter 
(Albright,  1980).  Extreme  surface  winds  of  20  m/s  are  observed.  (Brower 
et  al.,  1977).  Mean  of  vector  winds  over  days  105-180  is  71  cm/s  at 
065°. 

C)  Ice  Motions  AIDJEX  1972,  1975-6 

Mean  Drift  (summer)  mean  drift  is  ~  2.4  cm/s  (Hunkins,  1974)  in  the 
same  direction  as  the  geostrophic  current. 
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Fluctuating  ice  speeds  are  ~  10  cm/s  with  extremes  to  36  cm/s  in 
summer  and  less  in  winter  ~  1  cm/s  with  extremes  to  ~  20  cm/s  in  the 
CAB.  One  station  moved  140  cm/s  through  highly  fractured  ice.  (Thorn¬ 
dike  and  Colony,  1980). 

Mean  drift  in  winter  is  about  the  same  -  2  cm/day.  Fluctuations  in 
winter  are  less  -  1  cm/s  with  extremes  to  20  cm/s  (Thorndike  and  Colony, 
1980) . 

Near  shore  ice  is  essentially  motionless  all  winter. 

D)  Ice  cover  (Weeks,  1976) 


Mean 

Range 

Offshore 

Summer 

78% 

8  -  100% 

Winter 

99% 

70  -  100% 

CAB 

Summer 

92% 

30  -  100% 

Winter 

99% 

98  -  100% 

E)  Ice  Thickness  Distribution  (topography) 

Yearly  average  from  Blidberg  (1979)  of  all  AIDJEX  data  gives  an 
average  thickness  h  =  426  cm  . 

In  winter  Swithinbank’ s  (1971)  and  Le  Schack's  (1972)  thickness 
distributions  give  an  average  of  h  =  628  cm. 

In  August  1963  Thorndike  et  al. ,  (1974)  estimate  an  average 
thickness  distribution  of  h  =  330  cm. 

•  • 

F)  Dissipation  of  energy  by  ice  deformation  in  winter  =  +  r 

also  depends  on  strength  (Table  1)  and  location  since  deformation  is 
localized  0-20  m  in  middle  of  park  to  20  -  100  m-  — —  in  boundary 

shear  zone  (PritcharS,  1981).  m 

In  summer  more  than  10%  (usually  10  -  30  %)  of  the  pack  is  open  water, 
though  strains  may  be  large  the  stresses  are  small.  Energy  dissipation 
may  occur  through  uncorrelated  high  frequency  motions  causing  collisions 
between  ice  sheets.  Motions  between  the  AIDJEX  stations  BB  and  SB  are 
essentially  uncorrelated  at  frequencies  above  1/8  per  hour  and  have  mean 
velocities  of  -  .5  cm/s  (Colony,  1980  unpub.)  from  spectral  density 
plot.  A  flow  grinding  against  another  could  cause  an  appreciable  release 
of  energy  locally.  The  strength  of  pack  ice  in  summer  has  been  shown  to 
be  negligible,  however,  so  that  the  energy  released  overall  by  this 
mechanism  should  be  negligible  as  well. 
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Table  1.  Crushing  and  Buckling  Strengths  from  Rothrock  (1980) 


Ice 

ickness 

Crushing 

Stress 

Buckling 

Stress 

Ridging 

Stress 

n 

105  N/m 

10^  N/m 

105  N/m 

0.10 

0.3 

0.2 

0.0017 

0.50 

4.5 

2.7 

0.042 

1.00 

15.0 

8.7 

0.17 

3.00 

60.0 

50.0 

1.5 

G)  Drag  Coefficients 

Ice-water  drag  coefficient  for  lm  reference  velocity 


C  =  4.14  x  10  (Langleben,  1980) 


Stress  is  given  by 


o  C,  u-. 
w  1  1 


where  u^  is  current  velocity  at  1  m  below  ice,  Pw  =  density  of  water. 
Air-ice  drag  coefficient  for  10m  wind  velocity  C  n  =  1:58  x  10  3 


(Banke  et  al.,  1980)  for  pack  ice  free  of  major  ridges.  Leavit  (1980) 
gives  C^q  ~  2  x  10  3  . 

Air-sea  drag  coefficient  depends  on  sea  state.  Kraus  (1972)  gives 
C1Q  =  1.3  x  10"3. 
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